DISTRIBUTION OF APPROXIMANTS AND GEODESIC FLOWS 
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Abstract. We give a new proof of Moeckel's result that for any finite index subgroup of 
the modular group, almost every real number has its regular continued fraction approximants 
^jf^ equidistributed into the cusps of the subgroup according to the weighted cusp widths. Our 

■ proof uses a skew product over a cross-section for the geodesic flow on the modular surface. 

Our techniques show that the same result holds true for approximants found by Nakada's 
Q-continued fractions, and also that the analogous result holds for approximants that are 
algebraic numbers given by any of Rosen's A-continued fractions, related to the infinite family 
of Hecke triangle Fuchsian groups. 



1. Introduction 

For each nonzero real number, the regular continued fraction algorithm produces a sequence 
of rational numbers that are best approximations in the appropriate sense. It is natural to 
ask about the number theoretic properties of the numerators and denominators in each such 
sequence. In particular, as we learned from Sheingorn [She93j . going back to the 1940s there has 
been a study of the distribution into congruence classes of these numerators and denominators. 
The situation was beautifully resolved in 1982 by Moeckel |Moe82| . who used the celebrated 
connection between regular continued fractions and geodesies on the hyperbolic modular surface. 

A fascinating aspect of the problem is that it does not directly involve observables of the 
dynamical system defined by the Gauss map, that is by the interval map related to the regular 
continued fractions. Furthermore, the vocabulary of the result indicates some of its depth: (1) 
any finite index subgroup H of SL(2, Z) uniformizes the hyperbolic surface H\h\ where H is the 
Poincarc upper half-plane and H denotes the image of H in T := PSL(2, Z) = SL(2, Z)/ ± I; 
(2) a cusp k of H is an _ff-orbit in Q U {oo}; (3) the width of a cusp is the index w(k) of the 
5-stabilizer of a point in this i/-orbit inside the point's T-stabilizer. 

Combined with a later result of Nakanishi |Nak89j . the Moeckel result, that we reprove here 
— and then generalize — is the following. 

Theorem 1 (Moeckel, Nakanishi). Suppose that H is a finite index subgroup H C SL(2, Z). 
Then, for almost every real x, the regular continued fraction approximants of x are distributed 
in the cusps of H according to the relative cusp widths. That is, for each cusp k and almost all 

x, 

#{°<n<N M „e«} = 

N— >-oo iV 

Since for each positive integer m, the principal congruence subgroup r(m) defined as the 
kernel of the homomorphism to the projective matrix group over the finite ring Z/mZ is a normal 
subgroup, the result states that for almost all x, the approximants Pk/lk are equidistributed in 
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the cusps of r(m). In particular, when m = 2 the result is that for almost all x, its approximates 
have limiting frequency of 1/3 for each of the type odd/odd, odd/even and even/odd. 

Mocckel relied on the ergodicity of the geodesic flow on the unit tangent bundle of any finite 
cover of the modular surface. A few years thereafter, Jagcr and Liardet |JL88| gave a new proof 
(in the main setting of the principal congruence subgroups) that was in a sense purely ergodic 
theoretic; in particular, they used a skew product over the Gauss map to avoid the use of the 
geodesic flow. The logical structure of both proofs is of two steps: an ergodicity argument 
and a counting argument. Whereas Moeckel's counting argument seems slightly artificial - 
indeed, it forced a restriction to a particular subset of the finite covers of the modular surface, a 
restriction later removed by Nakanishi but only by the addition of a second counting argument 
in the complementary case — Jager and Liardet's counting argument is straightforward. On 
the other hand, the ergodicity argument in Moeckel's approach is in some sense immediate, but 
Jager and Liardet must labor in the corresponding step. 

We give a combined proof that unites the strengths of the two approaches. By using a skew 
product based over a cross-section of the geodesic flow on the unit tangent bundle of the modular 
surface, ergodicity is immediate and counting is straightforward. Furthermore, our proof needs 
no adjustment to treat the setting of any finite cover of the modular surface. The cross-section 
of the geodesic flow on the unit tangent bundle of the modular surface that we use is that found 
by Arnoux |Arn94j : this has the Gauss map as a factor transformation. To be more exact, the 
first return map to Arnoux's cross-section gives a double covering of the natural extension of 
the Gauss map. Our approach is quite flexible, as we show by not only recovering the Moeckel- 
Nakanishi results for the approximants generated by the regular continued fractions, but also 
extending these results to Nakada's a-continued fractions |Nak81| , and to Rosen's A m -continued 
fractions |Ros54j . These latter are related not to merely the modular group, but to each of an 
infinite family of triangle Fuchsian groups, known as the Hecke groups; the Rosen fractions give 
approximation by quotients of certain algebraic integers. 

The study of the intertwined nature of continued fractions and the geodesic flow on the mod- 
ular surface has a rich history. One of the first significant steps was in 1924, when E. Artin 
jArt24j gave a coding of geodesies in terms of the regular continued fraction expansions of the 
real endpoints of their lifts. A decade later, Hcdlund Hcd35 used this to show the ergodicity 
of the geodesic flow on the modular surface. A few years thereafter, E. Hopf proved ergodic- 
ity for the geodesic flow on any hyperbolic surface of finite volume, see his reprisal |Hop71 . 
Series |Ser85j (inspired by Moeckel's result) and Adler and Flatto, see especially their |AF91j . 
give cross-sections for the geodesic flow on the modular surface such that the regular continued 
fraction map is given as a factor. Quite recently, D. Mayer and co-authors [MS08] . |MM10j have 
given explicit cross-sections for geodesic flow on the surfaces uniformized by the Hecke triangle 
groups such that a certain variant of the Rosen continued fractions |Ros54] occurs as a factor. 
(The second-named author first learned of such possibilities from A. Haas, see the related treat- 
ment of another variant of the Rosen fractions in |GH961 .) For much of this history, further 
motivation, and also the work of S. Katok and co-authors using various means to code geodesies 
on the modular surface, see |KU07) . As stated above, for the regular continued fractions we 
use a cross-section given by Arnoux |Arn94j ; for the other continued fractions we study, we use 
cross-sections found by Arnoux and Schmidt jASj with a method directly inspired by |Arn94] . 
(Arnoux's approach to constructing cross-sections is more algebraic in flavor than the others we 
have mentioned.) 
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Our approach does have the disadvantage of presenting consecutive approximants as first or 
second columns alternately (see Lemma [5]), but this is natural when directly following the Gauss 
map. We address this by twisting from a conjugate of the subgroup of interest to the subgroup 
itself, see the end of the proof of Corollary [2j 

After this work was completed, we learned from Pollicott's recent survey article [Polll] that 
C.-H. Chang and D. Mayer |CM001 had previously used the skew product approach to lift a 
cross-section for geodesic flow on the unit tangent bundle of the modular surface to give cross- 
sections for finite index covers; this allows them to study dynamical zeta functions for these 
covers. 



2. Background 

2.1. Regular continued fractions. The regular continued fractions map, or the Gauss interval 
map, on D := [0, 1) is 

for x ^ ; T(0) := . 



x 



For put d(x) := [l/x\ , with d(0) = oo, and set d n = d n (x) := d(T n 1 (x)). This yields 

the regular continued fraction expansion of x G R : 

x = d H =: [do; di,d 2 , . . .] , 

d 2 H 

where do G 2 is such that x — do G D. (Standard convergence arguments justify equality of x 
and its expansion.) The approximants p n /q n of x G D are given by 



and 
(2) 





Po\ 







\q-i 


Soy 


-c. 





Pn-1 Pn\ = ( 1/0 1 \ 1 1 

q n -i q n J V 1 ^1/ V 1 ^2/ V 1 



for ?i > 1. Of course, this is equivalent to p n /q n = [di,d2, ■ ■ • ,d n ]. From the definition it is 
immediate that \p n -iq n — q n -iPn\ = 1 and in particular that the rational numbers p n /q n are 
all in reduced form. 

The standard number theoretic planar map associated to continued fractions is 

T(x, y) := (t(x), 1 ) , (x, y )en 
V d{x) + yj 

with O = x [0, 1]. The measure ijl on R 2 given by 

(3) dfx = (1 + xy)~ 2 dxdy 

is invariant for T. Indeed, Nakada, Ito, Tanaka |NIT77j showed that this two dimensional 
system gives a natural extension for the map T. In particular, the marginal measure obtained 
by integrating fi over the fibers {x} x {y \ (x, y) G fi} gives the invariant measure v with 
dv = (l + a;) -1 ^. The probability measure on II obtained by division by log 2 is known as Gauss 
measure. 
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2.2. Geodesic flow. Some of the following can be found in }Man91| . 

Using the Mobius action of SL,2(R) on the Poincare upper-half plane DH, one can identify a 
matrix with the image of z = i under the matrix. This results in SL,2(R)/S02(R) being identified 
with the Poincare upper half-plane H. Similarly, G := PSL2(R) = SL2(R)/ ±1 can be identified 
with the unit tangent bundle of DH. Recall that a Fuchsian group T is any discrete subgroup of 
G, cquivalcntly T acts properly discontinuously on H and the quotient is a hyperbolic surface 
(to be precise, possibly an orbifold). The unit tangent bundle of this quotient can be identified 
with T\G. 

The geodesic flow in our setting acts on the surface's unit tangent bundle: Given a time t and 
a unit tangent vector v, since the unit tangent vector uniquely determines a geodesic passing 
through the vector's base point, we can follow that geodesic for arclength t in the direction of 
v; the unit vector that is tangent to the geodesic at the end point of the geodesic arc is the 
image, gt{v), under the geodesic flow. The hyperbolic metric on DH corresponds to an element of 
arclength satisfying ds 2 = (dx 2 + dy 2 )/y 2 with coordinates z = x + iy. In particular, for t > 0, 
the points z = i and w = e t i are of distance t apart. Since G acts on the left by isometrics 
on H, the geodesic flow on its unit tangent bundle is given by sending A G G to AE t , where 
fe*/ 2 \ 

Ef = I _ t j 2 I . Similarly, on Y\G one sends the class represented by A to that represented 

by AEt ; we sometimes denote this by g t ([A]r) = [AE t ]r- 

There is a natural measure on the unit tangent bundle T DH: Liouville measure is given as 
the product of the hyperbolic area measure on H with the length measure on the circle of unit 
vectors at any point. This measure is (left- and right-) SL,2(R)-invariant, and thus gives Haar 
measure on G. In particular, this measure is invariant for the geodesic flow. When discussing 
measure theory, we always consider the Borel cr-algebra or its ^-completion. 

Let (X, n) be a measure space and a measure preserving flow on X, that is $ : X x R — >• 
X is a measurable function such that for $>t(%) = ^(x,t), $ s +f = $ s o $ t . Then E C X is a 
measurable cross-section for the flow $ t if: (1) the flow orbit of almost every point meets E; (2) 
for almost every x £ X the set of times t such that $t(x) G E is a discrete subset of R; (3) for 
every r > 0, the flow box ^4[o,r] = {®t(A) \ A G E, t <G [0, r] } is /i-measurable. This last defines 
the induced a- algebra on E. 

The return-time junction r = is r{x) = inf{< > : &t(x) G E} and the return map 
R : E — > E is defined by R(x) = ^ r ^(x). The induced measure fis on E is defined from flow 
boxes: one sets fJ,z(A) = ^fj,(A[ ^) for (all) < r < inf xe A{r(x)}. The flow (X 7 3§,[i) is 
then naturally isomorphic to the special flow built over the return map (E,/i£,i?) with return 
time r, as follows: define E ~ {(x,t) : < t < r(x)}/ ~, where ~ is the equivalence relation 
(x,r(x)) ~ (R(x),0), equipped with the measure /I, the product of /Lt£ on E with Lebesgue 
measure on R; the flow Q R on E is defined by $f : [x, s) n- (x, s + 1); this preserves ji. 

Convention: In all that follows, we will write cross-section to denote measurable cross-section. 

A flow $ t is ergodic if for any invariant set either it or its complement is of measure zero. 
Of fundamental importance for us is Hopf's result that the geodesic flow is ergodic on the unit 
tangent bundle of any finite volume hyperbolic surface, see his reprisal in [HopTI]. A flow is 
recurrent if the 4>-orbit of almost every point meets any positive measure set infinitely often. By 
the Poincare Recurrence Theorem, an ergodic flow on a finite measure space is recurrent. Given 
a cross-section, a return-time transformation (E, fj,s, i?s) is ergodic and recurrent if and 
only if the flow is. A semiconjugacy, homomorphism or factor map of two measure-preserving 
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flows or transformations is an a.s. onto map which preserves the measures and the dynamics; a 
conjugacy or isomorphism is in addition invertible, with the inverse map also a homomorphism. 
The following summarizes some of the above from a perspective of central importance to us. 

Lemma 1. Let G denote the group PSL(2, R). Suppose that T is a Fuchsian group of cofinite 
volume and that s/ C G projects injectively to give a cross-section E for the geodesic flow on 
the surface uniformized by T. Then for almost every A € s/ there exists a unique pair of first 
return time t = t^ and AI G T such that MAE t € s/ . 

We shall also need the following about lifts of cross-sections: 

Lemma 2. Let (X,£§,ji,(fr t ) and (X, S3, /j,, $j) be two measure-preserving flows such that the 
second is a factor of the first, i.e. there exists tt : X — > X measure-preserving and onto such 
that 7ro $ t = $ t o 7T. Let (E, R) be a cross-section and return map for <& t . Then E = 7r (E) 
is a cross-section for <J> i; with return function r = ro7r, return map R = ^^(x) and measure 
Jl on E the induced measure; the two return transformations are semiconjugate by way of the 
restriction of n to E . 

The proof follows from the definitions. We call (E, R) the lifted cross-section. 



2.3. Arnoux cross-sections. Inspired by Veech's "zippered rectangles" |Vee84| as a means to 
study the geodesic flow on Tcichmuller spaces, Arnoux |Arn94| gave a coding of geodesic flow on 
unit tangent bundle of the modular surface PSL(2, Z)\G by way of "boxes" described in terms 
of two types of subsets of SL2(R). Given (x,y) £ R 2 , let 



A- 1 (x,y)=(l x I and 



(4) 



Let 



A+i{x,y) = 



xy 

x 1 — xy 

-1 y 



(5) Z:R 2 \{(x,y)\y=-l/x}^R 2 

(x,y) i ^ (x,y/(l + xy)) 

and define the following subsets of SL2(K) (recall that f2 is the domain of the planar natural 
extension of the regular continued fractions map): 

A-i = {A-!(x,y) | (x,y) e ^(ft)} and A +1 = { A +X {x,y) \ (x,y) g ^(f!)} , 

and finally let s/ = A-i U A+i, considered as elements of G. 
Arnoux Ar n94] shows the following. 

Theorem 2. [Arnoux] The sets/ C G projects injectively to give a cross-section for the geodesic 
flow on the unit tangent bundle of the modular surface, thus on T\G = PSL(2, Z)\PSL(2, R). 
Furthermore, the first return map to this cross-section by the geodesic flow is given by the 
projection of 



A a {x,y) i ^ MA a (x,y)E t , 



(i 
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G 



M 




[M] r 



H\G 



[M] H 



Figure 1. Basic notation of covers: sets and element- wise maps. 



with t = — 21ogx and 



M = < 




1 



1 \l/x\ 



\l/x\ 1 



1 



if a = -1; 



ifa = +l. 



Moreover, the dynamical system defined by this first return map has the Gauss map as a factor 
via the projection to the first coordinate. 

2.4. Cusps, cosets and covers. An clement of G = PSL(2, R) is parabolic if it fixes a unique 



|trA| = \a + d\ being equal to 2. If T is a Fuchsian subgroup of G, then the T-orbit of any of its 
parabolic fixed points consists solely of parabolic fixed points for T; we call such an orbit a cusp 
of T. It is also a standard fact that the T-stabilizer of any such parabolic fixed point is infinite 
cyclic. By conjugating as necessary, we can assume that the parabolic fixed point is at infinity, 
and thus the stabilizer is a cyclic group of translations; the cusp of the group corresponds to 
a finite area end of the hyperbolic surface r\fH, traditionally also called a cusp — with our 
normalization, this geometric cusp is the projection of a hyperbolic triangle with one vertex at 
infinity, where the stabilizer acts so as to identify the two sides of the triangle. 

If H C r is any finite index subgroup, then any parabolic fixed point of V is also a parabolic 
fixed point of H. In particular, each cusp of T is partitioned into cusps of H. The width (with 
respect to T) of a cusp k of H is the index w(k) of the iJ-stabilizer of any representative for 
the cusp in the r-stabilizer of this representative parabolic fixed point. The geometric picture 
is that one can glue together w(k) copies of the triangle in H giving the cusp of the surface 
uniformized by T to construct the triangle giving the chosen cusp for the surface uniformized 



We denote a right ii-coset in T as either H-f or 7 mod H , depending upon ease of typography. 

Now, given a parabolic fixed point p of T and a coset H-f, since 7 • p is a parabolic fixed 
point (also) for H, it lies in some cusp k of H. But, then so does its whole f/-orbit; that is H-f 
sends p into k. There is thus an equivalence relation on the set of cosets H\T, characterized by 
sending p to the same cusp. Suppose that T has exactly one cusp. Since T acts transitively on 
its single cusp (viewed as a T-orbit), and T is the union of the right cosets of H in T, we see 
that every TJ-cusp k has a nonempty corresponding equivalence class. Let w'(k) be the number 
of elements in this class. 



element 




by H. 
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Choose a fundamental domain, for T acting on OH which reaches the boundary at p. Wc 
can also choose SeT generating the T-stabilizer of p. Upon choosing coset representatives ji , 
a fundamental domain for H is given by taking the union of these ji applied to & ' . Each ji in 
the equivalence class of k conjugates the generator S to a generator of the T-stabilizer of their 
common image of p. Indeed, the iT-stabilizcr of this image point is generated by the product 
over the equivalence class of the 7iS'7j~ • Therefore, we have the equality w'(k) = w(n). 

We summarize the above in the following. 

Lemma 3. Suppose that T is a Fuchsian group with a single cusp and that H is a finite index 
subgroup ofT. Then upon choice of a parabolic fixed point ofT, there is an equivalence relation 
on the right cosets H\T, corresponding to cusps of H . The size of an equivalence class is the 
width of the corresponding cusp. 

We define a notation for the equivalence classes on the right cosets. 

Definition 1. Given T and H as in Lemma El and a parabolic point p for T, let 

denote the equivalence class of the right coset Hj under the relation defined by equality of image 
of p, and 

denote the partition block of H\T corresponding to n under this equivalence relation. We refer 
to the relation itself as the -p relation on H\T. 

Note that we will occasionally slur over distinctions and refer to the various equivalence 
classes of H\T defined by a -p relation as the cusps of H . 

We will have need to switch roles between a group H and its conjugate by 

(6) fjso. 

Note that = u in G. 

Lemma 4. Suppose that T is a Fuchsian group with a single cusp, with iGT, and that H is 
a finite index subgroup ofT. //{7i}i<i<n is a full set of coset representatives for lHi\T, then 
the L~fi form a full set of coset representatives of H\T. Furthermore, the map p i— > i ■ p places 
the iHi-cusps in 1-to-l correspondence with the H-cusps. This correspondence preserves cusp 
widths. 

Proof. If r is partitioned by the iHu^i, then since t € V we have that r = iT is partitioned by 
the Hc^i. Thus the 17^ do give a full set of iJ-cosets. 

If p is a parabolic fixed point for V, then left multiplication by l sends lHl ■ p to H(l ■ p). If 
iHi^i ■ p = iHb'yj ■ p, then by left multiplication by l^ 1 — t, we immediately have H{b r fi) ■ p = 
H(L'jj) ■ p. Thus, the equivalence relation defined by p gives the same cusp widths for H as for 
lHl. 

□ 

For H a subgroup of G, we also use the notation for the associated projections shown in 
Figured] — in particular, for ease of reading, we often use an expression of the form [M]r to 
denote Trr(M). 
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3. Lifting cross-section as skew product 

Given a cross-section for the geodesic flow on a hyperbolic surface and a finite degree cover 
of the surface, from the general facts of Lemma [5] the cross-section can be lifted to give a cross- 
section for the geodesic flow on the cover; as we see here, in this case the resulting system can 
be nicely expressed as a skew product over the original cross-section. 

Theorem 3. Suppose that T is a Fuchsian group of cofinite volume and that the projection 
map 7tt injectively maps a subset si C G to a cross-section E for the geodesic flow on T\G = 
T 1 (L\IH). Let H C L be any finite index subgroup ofT, and H\T the set of right cosets of H in 
T. Then the skew product transformation 

S:#/X H\T si X H\T defined by 

(A, 7 modi?) i— > ( AdAEt, ^M" 1 mod H ) , 

where M and E t are as in Lemma [7J defines a dynamical system naturally isomorphic to the 
first return map on the canonical lifted cross-section of the geodesic flow on H\G = T 1 (lf\IH). 
The measures involved are these: on T\G and H\G we have Liouville measure, and their cross- 
sections E, E carry the natural induced measures; the measure on si is passed over via the 
bijection to £ and on the finite set H\T we take counting measure, with the skew product carrying 
the product of these. 

Sec Figure[2]for diagrams related to this proof. We use juxtaposition of sets of group elements 
to denote the set of all products of respective elements of these sets. 

Proof. Choose a set of right coset representatives ^ C F for H\T . 

Consider Ao £ si; its projection [^4o]r € irr(si) = E is sent to some [Ai]r by the first return 
map R on E. We claim that R is isomorphic to the map on si defined by Ao i-> MqAqE^ . To 
check this, by Lemma Q] the geodesic flow on T\G takes [j4o]r to [Ajjr, with A\ = MqAqEq , 
where as usual Mq G T and Eq is the diagonal matrix realizing the flow for a time of to . 
Next we consider 70 G ^ and 70^0 & \ we claim that the return map R for the lifted 
cross-section E is isomorphic to the map on ^si defined by 7oj4o ^ 71^1 where 71 = 7oM _1 . 
But since F uniquely decomposes as H'tf , there exists a unique pair ft e J?,7i € ^ such that 
/i _1 7i = 7oAf ~ 1 G r . Equivalently, hjo = 71 Mo and thus, 

hj A E = -/xMoAqEq = 71A1 . 

That is, [70^0 ]h returns to HH^si ) at time to . This is, moreover, the first return to this set: 
suppose that there are h' G H, 7 G 'if, A G si and t > such that h'joAoE t — -fA . The element 
in F given by M ~ 7~ 1 /i / 7o gives MAoE t = A and thus t is also a $r-return time for L4o]r , 
showing that t > to ■ Therefore, to is indeed the first return time for [70^-0] h- 

Since the map sending (A, Hj) to [7A]// is clearly surjective from si x H\T to tth ( ^ si ) , 
it remains to show it is also injective. Now suppose that A, A' G si and 7, 7' G T such that 
YiA] H = [l'A'] H - It follows that also [jA] r = [j'A'] T and hence ir r (A) = nr(A') ; by the 
injectivity of ir-p on si, we have A = A' . But, [7^]^ = [7'^']^ implies the existence of h G H 
such that hjA = j'A'; since A = A' we find that 7 and 7' are in the same right i7-coset of T. 

We have verified that our map is a bijection, which conjugates the skew product transfor- 
mation to the return map of E. The measures on the cross-sections correspond as well, since 
the Liouville measure on H\G may be viewed as the restriction, or projection, of Haar measure 
on G, while the induced measure on the cross-section E is the product of Haar measure on G 
restricted to si times counting measure on H\T (since a flow box for the lifted cross-section is 
a finite disjoint union of copies of the corresponding flow box for the section). 
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£/ x H\T x H\T (A, jmodH) h-^ (M AE t , -jM' 1 modiJ) 




[lAE t ] H 



AE, 



t jr 



Figure 2. Lifting the cross-section £ = 7rr( ). Diagrams of sets and element- 
wise maps. 



Finally, since the skew product is independent of choice of cosct representatives, we find that 
the isomorphism is natural. This completes the proof. □ 

Since invariant subsets of a cross-section and of a flow correspond, and since we know from 
Hopf's theorem that the geodesic flow on a finite volume hyperbolic surface is ergodic, the 
following is evident. 

Corollary 1. The skew product transformation S is ergodic. 

We note the following connection between flow cross-sections and fundamental domains. 

Proposition 1. If a G projects injectively to give a cross-section with return time r for the 
geodesic flow on the surface uniformized by T, then a (measurable) fundamental domain for the 
action ofT on G is given by the set 



{AE t | A e #f,0 < t < r} 



4. Distribution of regular continued fraction approximants into cusps 

Now we let T = PSL(2, Z) be the modular group, and H be any of its finite index subgroups, 
and begin with the cross-section stf for the geodesic flow on T\G given by Arnoux, see Theorem[2] 
above. The first return map to this cross-section maps elements of the (a = — l)-component to 
the (tr = +l)-componcnt and vice versa. This aspect of the map affects the presentation of the 
approximants Pk/qk- 

Lemma 5. Let &/ be the cross-section of the geodesic flow on the unit tangent bundle of the 
modular surface given in Theorem^ For A_i(x,y) £ , the second component of the k th 
composite of the skew product transformation S ( defined in Theorem 0) with itself applied to 
[A, I modTT) is 



proj 2 (S fc (A,/modtf)) = 



qk -Qk-i 
-Pk Pk-i , 



qk-x 



-qk 



-Pk-i Pk 



mod H if k is even; 



mod H if k is odd . 
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Proof. Since the initial value of a = — 1 and this value alternates in sign with each application 
of S, 

proj 2 (5' c (A,/ modff) ) = M^M^ 1 ■ ■ ■ M fc _1 mod H 

Upon referring to Equation ([2]), an elementary proof by induction verifies the statement. □ 

We need to address both the fact that the matrices above are of two types, and that the 
transpose of {pk,qk) is not a column of cither of the two matrices displayed above. We begin 
this with the following. (Recall that by Definition [TJ we can identify cusps with equivalence 
classes of H\T for the -p relation, where p is a parabolic fixed point of T.) In the following, we 
denote the indicator function for a set E by 1e, and we normalize the skew product measure to 
have total mass one. 

Lemma 6. Fix a cusp k of H , and let be the function 

* : srf X H\T -> R 
( A, 7 mod H ) H. (l^^xf^}^ + 1a +1 x{«}. )/2. 
Then, for almost every ( A, 7 mod H ) in stf X H\T, the limit 

1 N 

lim — V *o S k (A,-y mod H) 

fe=l 

is equal to the relative width of the cusp, w(k)/[T : H] . 

Proof. Since the function is a linear combination of indicator functions for sets of positive mea- 
sure, the crgodicity of S and the Birkhoff Ergodic Theorem implies the result. □ 

At this point, we can easily find that for almost every x, the series of values —qu/Pk derived 
from its regular continued fraction approximants has the appropriate frequency. To obtain the 
Mocckcl-Nakanishi result, we focus on the fact that i- {—qk/Pk) = Pk/qk- That is, we now finish 
our proof of that result. (When comparing with (p}, recall that we now insist on H denoting a 
projective group.) 

Corollary 2. For each cusp n of H and for almost every x £ D, we have 

V #{l<k<N\p k /q k eK} 

hm = w(k ) r : H\ . 

N— >oo IV 

Proof. The subset A-i x I mod H C srf x H\T has positive measure, and hence the limit of 
Lemma [6] holds for almost every element of it. But, this then implies that 

f[#7]oo if <r= -1; 
tf(4,(a:,2/),7niod.ff ) = < 

{[Hj}. ifa = +l. 

maps the corresponding sequence S k (A_i(x,y),I mod H) into the cusp k with limiting fre- 
quency equal to w(k)/[T : H]. From this it follows that for almost every 1 e 1, the sequence of 
—qk/Pk has this limiting frequency. 

We now use the above with the role of H replaced by lHl — we have that —qk/pk has this 
limiting frequency in any cusp of lHl. By Lemma [U we find that the sequence t • {—qk/Pk) 
has that same limiting frequency for corresponding cusps of H. By this same Lemma, the 
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cusp widths are the same for the two groups. But, conjugate subgroups have the same index. 
Therefore, the result holds. □ 

5. a-CONTINUED FRACTIONS 

Nakada [Nak81| defined the following family of interval maps. For a £ (0,1], we let D Q :— 
[a — 1, a) and define the map T a : D Q — > D Q by 

T a (x) 



1 






1 


+ 1 - a 










X 






x 





for x ^ 0; T Q (0) := . 



For x £ D Q , put 



e(x) := 



1 if x > , 
-1 ifa;<0, 



and d a (x) := 





1 


+ 1 - a 








X 





with d a (0) = oo. 

Furthermore, for n > 1, put 



s a ,n{x) ■= <T^-\x)) and d n = d a>n {x) := d a (T^-\x)). 



This yields the a- continued fraction expansion of x £ 

£ i 



?2 



di- 

«2i 

where do £ Z is such that a; — do £ Da- (Standard convergence arguments again justify equality 
of x and its expansion.) These include the regular continued fractions, given by a — 1 and the 
nearest integer continued fractions, given by a = 1/2. 

We define the a-approximants to x so that the recurrence relations 

p-i = 1; po = 0; p n = d n p n -\ + e n p n - 2 , n>l 

q-i = 0; g = 1; q n = d n q n -i + e n q n - 2 , n>l, 

hold. One easily checks that \p n -iq n — q r i-iPn| = 1- 

We show in this section that the distribution of the a-continucd fraction approximants pk / qk 
is the same as that for the regular continued fractions. 

Theorem 4. Suppose that H is a finite index subgroup H C PSL(2, Z). Then, for each a £ (0, 1] 
and for almost every real x, the a-continued fraction approximants of x are distributed in the 
cusps of H according to the relative cusp widths. 



The standard number theoretic planar map associated to these continued fractions is defined 



by 



T a (x,y) := T a (x) 



1 



(x,y) £ tt a 



d a (x) +e(x)y / 

see |KSS12] for a description of fl a C !„ x [0,1]. Here also, /i as given in is an invariant 
measure. 

Let \x a be the probability measure given by normalizing \i on fi Q , and v a the probability 
measure normalized from the marginal measure obtained by integrating /i Q over the fibers {x} x 
{y | (x, y) £ f2 Q }, 33 a the Borel er-algebra of D a , and £%' a the Borel tr-algebra of fl a . Kraaikamp- 
Schmidt-Steiner |KSS12| showed that (tt a , T a , 3$' al [la) is a natural extension of (B Q , T a , SS ai v a ). 

Analogously to the regular continued fraction case of a = 1, let 

Aa,-i = {A-x(x,y) | (x,y) £ Z~\n a )} and A +1 = { A +1 (x,y) \ (x,y) £ Z" 1 ^)}, 
where Z is given in J5|); and, let s/ a = A a -\ U A a ,+i, considered as elements of G. 
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Arnoux-Schmidt IASI show that 



Theorem 5. [Arnoux-Schmidt] The set si a C G projects injectively (up to measure zero) to a 
cross-section for the geodesic flow on the unit tangent bundle of the modular surface. Further- 
more, the first return map to this cross-section by the geodesic flow is given by the projection 
of 

A a (x,y)^ MA a (x,y)E t , 
with t = —2 log \x\ and M equaling 

l\ (1 d\ fd -1\ fl 
-1 d) ' ^0 V ' v 1 )\ d 1 

where d = d a (x), as (a,e(x)) = (— 1, — 1), (— 1, +1), (+1, — 1), (+1, +1), respectively. Moreover, 
the dynamical system defined by this first return map has the map T a as a factor. 

Note that when e(x) = +1, M is of the same shape as in the regular continued fraction case 
given by Arnoux, see Theorem[5] In all cases, si a has two components, indexed by a = ±1. In 
general, the map sends an element with e = +1 in the component indexed by a to the component 
indexed by —a, whereas elements with e = — 1 are mapped to within the same component. This 
is key to showing that the analog of Lemma [5] holds with only two types of matrices, instead of 
four as the previous result might suggest. 

Lemma 7. Let si a be the cross-section of the geodesic flow on the unit tangent bundle of 
the modular surface described in Theorem\^ Suppose that A-i(x,y) € stf a , then the second 
component of the k th composite of the skew product transformation S with itself applied to 
(A, I modH) is 



W 0) 2 (S k (A,ImodH)) 



qk -qk-i 
-Pk Pk-1 / 

q k -i -q k 

-Pk-i Pk , 



modH i/ proj 1 (5 fe (A, I modH)) e A a ,-i 



mod H otherwise. 



Note that the condition defining the two cases above becomes that of Lemma [5] upon restrict- 
ing to the case of all Ei(x) = +1. 

Proof. For A = A-\{x,y) £ A a ,-i, we have proj 1 (5 fc (^4, / modH)) € A a ,-i if and only if 
the number of £i(x) = +1 for 1 < i < k is even. The direct calculation verifying Lemma [5] 
shows that each application of S with £i(x) = +1 changes between the two shapes of matrices 
shown. Another direct calculation shows that each application of S with Si (x) = — 1 preserves 
the respective shapes. The result follows. □ 

It is now immediate that the analogous statements to Lemma [5] and Corollary [5] hold, and 
we indeed conclude that the limiting distribution into cusps of the a-approximants is the same 
as that for the regular continued fraction approximants. 

6. Rosen continued fractions 

Let A = A m = 2cos^-, then F m , the Hecke group of index m is the single-cusped Fuchsian 
group generated by the elements i given above in ([6]), and 



7m = 



1 A„ 
1 
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The group T3 is the modular group PSL(2, Z). For m ^ {3, 4, 6}, the group T m is non-arithmetic, 
see }Leu67] : and, thus in particular has no finite index subgroup that is PSL(2, R)-conjugate to 
a finite index subgroup of PSL(2, Z) . 

In 1954, D. Rosen |Ros54| defined an infinite family of continued fraction algorithms to aid 
in the study of the Hecke groups. The Rosen continued fractions and variants have been of 
recent interest, leading to results especially about their dynamical and arithmetical properties, 
see [BKSOOj , |NaklOj , }DKS09j ; as well on their applications to the study of geodesies on related 
hyperbolic surfaces, see |SS95j . [BS04] . |MS08I ; and to Teichmiiller geodesies arising from (Veech) 
translation surfaces, see [AHOOj . |SU10] . [5Ullj and |AS09j . Several basic questions remain open, 
including that of arithmetically characterizing the real numbers having a finite Rosen continued 
fraction expansion, see jLeu67j . |HMTY08] and |AS09] . 

We show in this section that the Rosen continued fraction approximants pu/qk are distributed 
analogously to those for the regular continued fractions. (We note that [NaklOj gives the fol- 
lowing result in the setting of the distribution into cusps of the principal congruence subgroups 
of the Hecke groups, using an approach directly related to Moeckel's.) 

Theorem 6. Fix an integer m > 3, and suppose that H is a finite index subgroup of the Hecke 
group of index m, H C r m . Then, for each for almost every real x, the Rosen X m -continued 
fraction approximants of x are distributed in the cusps of H according to the relative (to Y m ) 
cusp widths. 

Let D m = [—A/2, A/2) for m > 3. For a fixed integer m > 3, the Rosen continued fraction 
map is defined by 

fill " AL|^| + |J ^0; 

T m (x) = < 

[0 £ = 

for x € D m ; here and below, we omit the index "m" whenever it is clear from context. For n > 1, 
we define 

e n (x) = siT^x) and r n {x) = r(r^ 1 x) 

with 

e(y) = sgn(y) and r(y) 
Then, as Rosen showed in |Ros54| . the Rosen continued fraction expansion of x is given by 
[ei(x) : ri(x), e 2 (x) : r 2 (x), ... , e n (x) : r n (x), ...}:= ^-g- . 





1 


1 






+ 2 



r 2 A + 



As usual we define the approximants p n /q n of x € D m by 

( V-\ Pa\ _ A 



q-i qoj V° 1 

and 

r Pn-l Pn\ _ [0 £ l \ f° £ 2\ f0 £ n 



\ 1 \r\ ) V 1 Ar 2 / VI Ar„ 

for n > 1. From this definition it is immediate that \p n -iq n — q n -iPn\ = 1, and that the 
following well-known recurrence relations hold: 

P-i = 1; Po = 0; p n = Xr n p n -i + e n p n - 2 , n > 1 
q-i = 0; qo = 1; q n = Ar„g„_x + £ n g„_ 2 , n > 1. 
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Sketch of proof of Theorem® For each m, Burton-Kraaikamp-Schmidt BKSOO determined a 
planar natural extension on a region f2 m with the measure fi as above. Using this, Arnoux- 
Schmidt |ASj determine an stf m C G = PSL(2, R) such that £# m projects to a cross-section of 
the geodesic flow on T m \G. The first return map to this cross-section is a double cover of the 
natural extension. Furthermore, for each index m an analog of Theorem [5] holds as |AS| show, 
where again when Ei(x) = +1 there is a change of component, and when £i{x) = — 1 there is 
not; the matrices giving the return to the cross-section here are 

1 \ fl d\ m \ fd\ m -l\ / 1 N 



1 dX m J ' \0 1 J ' \ 1 J ' \dX m 1 

where d = d a (x), as (cr, e(x)) = (— f , — f ), (— f , +1), (+1, — f ), (+1, +1), respectively. 

Therefore, we easily find that the analog of Lemma [7] holds. Now, for each m we have both 
that 7 m fixes oo, and l 6 T m , this allows the proof of Corollary [2] to apply. From this, one easily 
completes the proof of Theorem [5] □ 



As an example of an implication of this, we consider the case of q = 5. Note that A5 is the 
golden ration, (1 + \/5)/2. 

Example 1. Almost every real number has Rosen As-continued fraction approximants with 
equal asymptotic frequency of the five types: odd/even, even/odd, 1/1, I/A5, I/XA5 + 1). 

The proof of this relies on checking the data related to the principal congruence subgroup of 
G5 for the ZfAs] ideal generated by 2 — this (normal) subgroup has five cusps, each of width 
two. We have chosen cusp representatives in the (75t)-orbit of 00. 

We briefly sketch a non-geometric manner to see that there should indeed be five cusps for this 
principal congruence subgroup. Note that Z[As] is the full ring of integers of its quotient field, 
and that the rational prime ideal (2) remains prime (is inert to this number field extension). 
Thus, each element of a pair (pk,Qk) reduces modulo (2) to some element of the quotient field 
Z[As]/(2), a finite field with four elements. However, each pair is relatively prime, thus there are 
obviously at most twelve combinations possible. Now, if either of Pk,qk is a multiple of 2, then 
multiplying each by an inverse modulo 2 of the other reduces the quotient to the form of one of 
the cusps odd/even, even/ odd. For the six remaining possibilities, one multiplies through by a 
multiplicative inverse of pk to see membership in one of the three remaining cusps. 

7. Further remarks 

We note that Theorem [3] remains true upon replacing the geodesic flow by the flow defined 
by any one-parameter subgroup of G; if this new flow is represented by the right multiplication 
by matrices of the form F t , then Lemma Q] and the theorem simply require the replacement of 
the E t by these F t . Of course the ergodicity of the skew product transformation announced in 
CorollaryQ] depends on the flow at hand. For the (stable and unstable) horocycle flows, results 
of Hedlund [Hed36j and later of Dani and Smillic |DS84j can be invoked, so ergodicity (and even 
essential unique ergodicity; that is, unique ergodicity after the suppression of weights on the 
collection of periodic orbits) of the skew product again holds. 
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